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Abstract. HE 0141-3932 (Zcm = 1.80) is a bright blue radio-quite quasar with an unusually weak Lya emission line. Large 
redshift differences (Az = 0.05) are observed between high ionization and low ionization emission lines. Absorption systems 
identified at Zabs = 1.78, 1.71, and 1.68 show mild oversolar metallicities (Z ~ 1 — 2Zq ) and can be attributed to the associated 
gas clouds ejected from the circumnuclear region. The joint analysis of the emission and absorption lines leads to the conclusion 
that this quasar is seen almost pole-on. Its apparent luminosity may be Doppler boosted by 10 times. The absorbing gas shows 
a high abundance of Fe, Mg and Al ([Fe, Mg, Al/C] ~ 0.15 ± 0.10) along with underabundance of N ([N/C] < -0.5). This 
abundance pattern is at variance with current chemical evolution models of QSOs predicting [N/C] > and [Fe/C] < at 
Z ~ Zq . 
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1. Introduction 

In the course of a high-resolution study of the Lya forest at 
intermediate redshifts (1.5 < z < 2) in bright quasars from 
the Hamburg/ESO Survey with the UV-visual echelle spectro- 
graph (UVES) at the VLT, the discovery spectrum of the QSO 
. HE 0141-3932 (z ~ 1.8, B = 16.2, Wisotzki et al. 2000) at- 
' traded our attention by two facts: it appeared to have no or only 
very weak Lya emission along with clearly recognizable Mg II, 
and it showed several absorption line systems at Zabs ranging 
from 1.78 to 1.68 with very complex and strong metal profiles 
which allow us to suggest that these systems may originate in 
the ejected gas. 

Since the identification spectrum of HE 0141-3932 was of 
only moderate quality, we took further low-resolution spectra 
with EFOSC 2 with the ESO 3.6 m telescope to improve the 
redshift measurement and we found a thkd peculiarity, namely 
large redshift differences (Az ~ 0.05) between high ionization 
and low ionization emission lines. 

Emission line spectra similar to that of HE 0141-3932 are 
observed in some high redshift BL Lacertae (BL Lac) objects 
(Urry & Padovani 1995, hereafter UP95). However, HE 0141- 
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3932 is a radio-quiet quasar, i.e. not a blazar A few other blue 
radio-quiet quasars are known with apparently missing or weak 
Lya emission but clearly present metal lines: PG 1407H-265 
[zem= 0.94, McDowell et al. (1995)], Tol 1037-2703 [Zem = 
2.20, Srianand & Petitjean, (2001)], PHL 1811 [Zem= 0.192, 
Leighly, Halpern, & Jenkins (2004)], and several cases from 
SDSS (Fanetal. 1999, 2003). The nature of their peculiar emis- 
sion line spectra is not clear 

In particular, Leighly et al. (2004) suggest a high accre- 
tion rate which powers the UV emission from an optically 
thick accretion disk, while suppressing the formation of a hot 
corona. However, recent radio observations of PG 1407H-265 
on the milliarcsecond scale revealed a relativistic jet of mod- 
erate power beamed toward us (Blundell, Beasley & Bicknell 
2003). This means that the quasar is seen almost pole-on and its 
emission line spectrum may be diluted by continuum radiation 
both from the jet and the accretion disk. 

In the present paper we analyze both the emission and ab- 
sorption line spectra of HE 0141-3932. We expect that this 
joint consideration may shed some light on the unusual proper- 
ties of this quasar. 

The paper is organized as follows: observations are de- 
scribed in Sect. 2, the emission lines are studied in Sect. 3, 
the analysis of the Zabs ~ ^om absorbers is given in Sect. 4, the 
results are discussed and summarized in Sect. 5. 
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Fig. 1. EFOSC 2 spectra of HE 0141-3932. The spectral resolutions are indicated in the panels. AH identified emission features 
are labeled. The parameters of the prominent emission lines are listed in Table 2 



2. Observations 

HE 0141-3932 was observed with the UVES at the 8 m ESO 
VLT on Paranal over 7 nights in July /August 2001. Eleven in- 
dividual exposures with integration times of 60 min were made 
using the dichroic mode in standard settings (Table 1). With a 
slit width of 1", a resolution of 41000 (7.3 km s~^ ) in the blue 
and 38000 (7.9 km s~^ ) in the red is achieved. 

The data reduction was performed at the Quality Control 
Group in Garching using the UVES pipeline Data Reduction 
Software (Ballester et al. 2000), the vacuum-barycentric cor- 
rected spectra were co-added. The resulting signal-to-noise ra- 
tio, S/N, is typically 75. 

Since both the Lya and the Mgii emission line ranges 
are covered only by UVES spectra, we performed an ab- 
solute flux calibration using an appropriate master response 
curve and the respective airmasses during the observa- 
tions. The procedure is described on the ESO web page 
www.eso.org/observing/dfo/quahty/UVES/qc/response.html. 

Further spectra were taken with EFOSC 2 at the ESO 3.6 m 
telescope on October 2, 2003, to improve our knowledge about 
the redshifts of the emission fines. Details are given in Table 
1 . Wavelength and flux calibration was performed according to 
standard procedures. The resulting spectra are shown in Fig. 1 . 
Notice that due to the rapidly decreasing sensitivity below 3500 
A, the EFOSC 2 spectra only partially cover the wavelength 
range expected for the Lya emission line. 



3. The emission iine spectrum 

The combination of EFOSC 2 spectra with roughly flux cali- 
brated UVES spectra allows us to estimate redshifts and equiv- 
alent widths of all emission lines between Lya and Mg II A2800 
A. Equivalent widths (EW) and redshifts were measured by 
fitting Gaussian profiles to the data. Table 2 presents the re- 
sults. The Lya EW is relatively uncertain since in the caU- 
brated UVES data the continuum is difficult to determine and 
the blue wing of the line is incomplete (see Fig. 1 and 2). 

The identification of the emission fine at ~ 3400 A is am- 
biguous: it could be either Lya at 2; = 1.80 or Nv at 2: = 1.75 
or blend of both lines. There are two arguments in favor of the 
N V identification: it has a width comparable to C iv and its red- 
shift - if identified as N v - is equal to that of the C iv, while in 
QSOs with significant redshift differences between high ion- 
ization lines (Civ) and low ionization lines (Mgll, H^), the 
Lya line typically is found at the redshift of the high ionization 
BLR emission fines (Gaskell, 1982; Espey et al., 1989). On the 
other hand, the emission line at 3400 A is quite strong for N V, 
although there are quasars known with N v strength compara- 
ble to that of C IV (Hah et al. 2004; Baldwin et al. 2003b). Even 
if this line would be entirely due to emission of neutral hydro- 
gen, we can conclude that Lya in HE 0141-3932 is unusually 
weak. 

The second anomaly of this quasar is the large velocity 
difference between low ionization lines (Mgii, Cm], Aim 
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Table 1. Log of spectroscopic observations of HE 0141-3932 



Instrument 


AA 


Exposure 


Date 


Resolution 




(A) 


time (h) 




(kms^^) 


UVES 


3053-3872 


7 




7.3 




4700-5800 


7 


July + 


7.3 




5800-6800 


7 


August 


7.8 


UVES 


3761^982 


4 


2001 


7.3 




6650-8550 


4 




7.8 




8600-10400 


4 




7.8 


EFOSC 2 


3270-5240 


0.25 


Oct. 2 


420 (6 A) 




3380-7520 


0.25 


2003 


700 (13 A) 



etc.) and high ionization lines (C IV, Si IV) of roughly 5000 
km s~^ ; it is among the most extreme cases like PG 1407+265 
(McDowell et al. 1995, AwMgii-civ ^ 5000 km s'^ ) and Q 
1309-056 (Espey et al. 1989, AuMgii-civ ^ 4000 km s'^ ). 

In Fig. 2 we display the emission line profiles on a velocity 
scale relative to 2; = 1.80. The redshift of the low ionization 
lines of z = 1.80 appears to be roughly the systemic redshift, 
since the UVES spectrum shows that the Lya forest starts at 
z = 1.8086. 

While the velocity shift in the emission lines is extremely 
large, how unusual are the emission line intensities ? This can 
be best seen by a comparison with the histogram of equivalent 
widths of the LBQS-quasars, a well selected sample of QSO 
given by Francis (1993). The Lya, C iv, C iii] and Mg 11 equiv- 
alent widths are among the bottom 2%, 1%, 13% and 1.5%, 
respectively, of the Francis (1993) distribution. At the same 
time, the Aim 1856/64 and Feiii 2075 lines appear unusually 
strong. In order to investigate the physical properties of gas 
which could produce such an unusual emission line spectrum 
we compared the observed equivalent widths of HE 0141-3932 
with the compilation of quasar BLR rest-wavelength emission 
line equivalent widths given as functions of column density, in- 
cident ionizing spectrum, and metal abundance of the emitting 
gas clouds (Korista et al. 1997). While Lya, Mgii and Cm] 
have roughly the same velocity (redshift) and could originate in 
the same volume, we were unable to find a parameter combina- 
tion which met the constraints on Lya and Mg 11 line strengths 
simultaneously - the theoretical Lya: is always much too strong 
for a parameter combination that matches the Mg 11 equivalent 
width. The C iv/Lya EW ratio can be met assuming high gas 
densities [log n(H) > 12], but the discrepant redshifts exclude 
formation in the same location. A mean BLR parameter set ap- 
pears to be unable to reproduce the measured EWs of emission 
Unes in HE 0141-3932. 

4. The absorption systems 

4. 1. Calculation procedure 

In order to estimate the physical parameters of the absorption 
systems we used the Monte Carlo Inversion (MCI) method. 
Detailed description of the MCI is given in Levshakov et al. 
(2000, 2002, 2003a). Here we outline briefly its basics needed 
to understand the results presented. 
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Fig. 2. Emission lines of HE 0141-3932 (spectra from UVES 

and EFOSC 2). v = km s~^ corresponds to z = 1.80. For 
orientation, alv = km s^^ and v = —5000 km s^^ dashed 
lines are plotted 



The MCI is based on the assumption that all Unes observed 
in the absorption system are formed in a continuous absorbing 
gas slab of thickness L where the gas density, nii{x), and ve- 
locity, v{x), fluctuate from point to point giving rise to complex 
profiles (here x is the space coordinate along the line of sight). 

We also assume that within the absorber the metal abun- 
dances are constant, the gas is optically thin for the ionizing 
UV radiation, and the gas is in thermal and ionization equilib- 
rium. The intensity and the spectral shape of the background 
ionizing radiation are treated as external parameters. 
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Table 2. Estimated equivalent widths, redshifs and relative velocities for emission lines of HE 0141-3932 



Ion 


Instrument 






V, km s ^ 


Comment 


Hi 1215 
and/or Nv 1239 


UVES 


16 ± 7 


1.75 ±0.01 


'-'-2140 

-5350 ± 1070 


asvmmetric 


Si IV 1398+Oiv] 1402 


EFOSC 2 


11 ± 1 


1.76 ± 0.01 


-4280 ± 1070 


— 


Civ 1549 


EFOSC 2 


Mil 


1.75 ± 0.01 


-5350 ± 1070 


asymmetric 


Aim 1858 + 

Cm] 1909 + Si III] 1892 


EFOSC 2 


13 ± 1 


1.79 ± 0.01 
1.79 ± 0.01 


-1070 ± 1070 
-1070 ± 1070 


blend, fractions of components unknown 


Feiii 2075 


EFOSC 2 


5± 1 


1.80 ± 0.01 


± 1070 




Mg II 2800 


UVES 


16 ± 1 


1.795 ± 0.005 


-535 ± 535 


blended with Fe ii emission 



Zero velocity corresponds to z = 1.80. Only statistical errors are indicated. 



The radial velocity v{x) and gas density nu{x) are con- 
sidered as two continuous random functions which are repre- 
sented by their sampled values at equally spaced intervals Ax. 
The computational procedure is based on adaptive simulated 
annealing. The fractional ionizations of different elements at 
each space coordinate x, Tion[U{x)], are computed with the 
photoionization code CLOUDY (Ferland 1997). 

The following physical parameters are directly estimated 
by the MCI procedure: the mean ionization parameter Uq, the 
total hydrogen column density A'h, the line-of-sight velocity 
dispersion, a^, and density dispersion, ay, of the bulk material 
[y = nB_{x)/no], and the chemical abundances of all ele- 
ments involved in the analysis. With these parameters we can 
further calculate the column densities for different species N^,, 
and the mean kinetic temperature Tkin- 

In general, the uncertainties on the fitting parameters Uq, 
Nn, (Ty, (7y, and Za are about 15%-20% (for data with S/N 
> 30) and the errors of the estimated column densities are less 
than 10%. However, in individual absorption systems the accu- 
racy of the recovered values can be lower for different reasons 
Uke partial blending of the line profiles, saturation, or absence 
of lines of subsequent ionic transitions. 

The MCI can be supplemented with an additional proce- 
dure aimed at restoring the shape of the background ionizing 
spectrum. A formal description of this procedure is given in 
the Appendix. Its accuracy depends significantly on the num- 
ber of unsaturated lines of the subsequent ionic transitions of 
different elements (e.g. Si Il/Si Ill/Si IV, C Il/C Ill/C IV) avail- 
able for the analysis. The absorber at Zabs = 1.7817 (described 
below) reveals enough such lines and we use them to estimate 
the spectral shape of the background UV radiation in the range 
E > 1 Ryd. 

All calculations are carried out with the laboratory wave- 
lengths and oscillator strengths taken from Morton (2003). 
Solar photospheric abundance for carbon is taken from AUende 
Prieto et al. (2002), for siUcon, nitrogen and iron - from 
Holweger (2001). 

4.2. Absorption system at Zabs = 1.7817 

Many unsaturated lines of different ionic transitions along with 
a single saturated hydrogen line Lya are detected in this system 



(Fig. 3) separated by 2 000 km/s from the central source. The 
blue wing of Lya is partially blended by the adjacent hydrogen 
line [iV(H l) ^ 10^'* cm^^ ] from an intervening absorber (only 
a weak C iv doublet is detected in this system). C iii A977 is 
beyond the wavelength coverage and Sim A1206 is blended, 
probably with Lya absorption from the ^abs = 1.7606. 

We started the analysis assuming standard ionizing back- 
grounds such as a power law oc j^^" (with different indexes 
a = 1.0 — 1.8), the mean intergalactic spectrum (at z = 1.8) 
of Haardt & Madau (1996), and the AGN-type spectrum of 
Mathews & Ferland (1987, hereafter MF). 

None of these spectra was consistent with the observed in- 
tensities of the C II, C IV, Si ii and Si iv Unes: all trials under- 
produced C n/C IV and overproduced Si ii/Si iv, but the MF 
spectrum provided the lowest x^- All spectra gave the mean 
ionization parameter in the range of —2.75 < log ?7o < — 2.5. 
In spite of the saturation, the neutral hydrogen column density 
can be estimated with a sufficiently high accuracy (~ 20%) 
since the velocity dispersion of gas is determined by numerous 
metal lines detected in this system (the procedure to restore 
partly blended profiles is described in Sect. 3 in Levshakov 
et al. 2003a). AH runs with different model spectra showed a 
rather high metallicity - sUghty above solar value, and an al- 
most solar ratio of Si/C. 

Taking into account this preliminary information we can 
assume that the UV spectrum to be found should maximize 
the product of ratios C Il/C IV and Si iv/Si ll in the range of 
Uq ~ 0.002 — 0.003 for solar metal content and solar element 
abundances. The spectral shape of the MF spectrum can be 
taken as a first approximation (for more details see Appendix). 

Applying the adjustment procedure as described in 
Appendix, we estimated a new shape of the ionizing spec- 
trum which ensured a self-consistent description of all lines 
observed in the ^abs = L7817 system. This spectrum is shown 
in Fig. 4 by the dashed line, whereas the initial MF spectrum is 
the soUd line. As seen in Fig. 4, the restored spectrum is softer 
in the range 2 Ryd < E < 6.5 Ryd but much harder above 
6.5 Ryd. It shows a break at the He ii ionization edge with the 
amplitude J^(A)/Jy(B) ~ 3 corresponding to the optical depth 
rc(He II) ~ 1. Physical parameters estimated with this modified 
ionizing background are listed in Table 3, Col. 2, and the corre- 
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Fig. 3. Hydrogen and metal absorption lines associated with the Zabs= 1.7817 system toward HE 0141-3932 (normalized inten- 
sities are shown by dots with 1 a error bars). The zero radial velocity is fixed at 2; = 1.7817. Smooth curves are the synthetic 
spectra convolved with the corresponding point-spread spectrograph function and computed with the physical parameters listed in 
Table 3, Col. 2. Bold horizontal lines mark pixels included in the optimization procedure. The synthetic profiles of unmarked ab- 
sorption features were calculated in a second round using the velocity v{x) and gas density nii{x) distributions already obtained 
in the optimization procedure. 



sponding synthetic profiles are plotted in Fig. 3 by the smooth 
curves. 

4.3. Absorption system atz^hs = 1-7103 

This absorption system spans the velocity range of 500 km s~^ 
and is located in velocity space at a distance of 9 600 km s~^ 
from the QSO. It reveals many lines of different ionic transi- 
tions (Fig. 5). The structure of the Lya, Civ, Siiv and Nv 
profiles indicates that the system may be divided into two sub- 
systems: one at -200 km s~^ < v < 180 km s~^ (A) and 
another at 180 km s"^ <v < 320 km s"^ (B). 



The subsystem B has an unsaturated hydrogen line and 
hence allows us to estimate the hydrogen column density with a 
sufficiently high accuracy. The physical parameters computed 
with the UV background deduced for the Zabs= 1-7817 sys- 
tem are given in Table 3, Col. 4. The obtained extremely high 
metallicity - almost 9Zq - is striking. Since silicon lines Si III 
A1206 and SilV A1393, 1402 are very weak and the element 
ratios ([Si/C], [N/C]) not known a priori, the mean ionization 
parameter Uq ~ 0.04 represents a lower limit consistent with 
the observed intensities of the C iv lines and the absence of the 
C II A1334 absorption. Formally this subsystem could be fit- 
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Fig. 4. Spectral shape of a typical AGN ionizing continuum 
from Mathews & Ferland (1987) shown by the solid line and 
its modification (dashed line) estimated for the 2:abs= 1.7817 
system. The spectra are normaUzed so that Jv{hy = 1 Ryd) = 
1. Factors fi are defined in the Appendix 

ted with any higher ionization parameter producing even higher 
metaUicity and higher [Si/C] and [C/N] ratios. 

To confirm the high metaUicity in this subsystem we re- 
peated calculations assuming other ionizing backgrounds - 
power laws with a ranging between 1.0 and 1.8 and the MF 
spectrum. These trials also delivered metaUicity of about 10 
solar. 

In principle, such high metal enrichment of the circumnu- 
clear gas is predicted in some models of chemical evolution of 
QSO/host galaxies (Hamann & Ferland 1999, hereafter HM). 
However, another explanation of our results is that the absorb- 
ing gas is not in equilibrium with the ionizing background, e.g., 
is stiU cooUng and recombining. In this case high metallicities 
can be caused by a longer recombination time of hydrogen as 
compared to C iv and N v (e.g., Osterbrock 1989). The ioniza- 
tion parameter U and, hence, the total hydrogen column density 
TVh oc 1/Thi(J7) is determined by the observed line intensi- 
ties of different ions. Since hydrogen is ionized higher than it 
would be in the ionization equiUbrium with the same param- 
eter U, the total hydrogen becomes underestimated leading to 
an artificially high metaUicity. We may expect that the adjacent 
subsystem A (centered at w = km s^^ ) would help to clarify 
the physical conditions in the .Zabs = 1-7103 absorber. 

The subsystem A shows complex profiles of low ionization 
lines C II, Mg II, Si II, and Fe II along with Al III, Si III, Si IV, 
N V and saturated C iv (see Fig. 5). The computations had been 
carried out with the ionizing background from the Zabs =1.7817 
system. The presence of the subsequent ionic transitions guar- 
antees the accuracy of the mean ionization parameter of ~ 10% 
(for a given ionizing background). 

The assumption of constant metalhcity inside the absorber 
turned out to be inconsistent with the observed blue wing of 
Lya and the absence of any metal absorption in the range —200 



km s~^ < V < —100 km s~^ . A single available hydrogen 
line does not allow us to conclude whether the blue wing is 
blended or there is indeed a gradient of metal content. The for- 
mer possibility looks more probable since there is a non-zero 
flux in the Lya profile at -110 km s~^ < v < -90 km s~^ 
with the mean intensity 0.018 ± 0.005 (marked by the arrow 
in Fig. 5). The constant metaUicity throughout this sub-system 
seems to be appropriate since very steep blue wings of several 
ions (Si III, Si IV,C IV) do not indicate a progressive dilution. 
Thus, the calculations were carried out with the assumption of 
constant metaUicity. The neutral hydrogen column density and 
the shape of the blue wing of the synthetic Lya shown in Fig. 5 
were calculated with the velocity and density fields estimated 
from the observed metal line profiles (for details see Sect. 4.3 
in Levshakov et al. 2003a). The obtained physical parameters 
are given in Table 3, Col. 3. We also tried other ionizing back- 
grounds (different power laws and the MF spectrum), but none 
of them was consistent with the observed relative intensities 
within the C Il/C IV and Si Il/Si Ill/Si IV profiles. 



In spite of all uncertainties intrinsic to this subsystem sev- 
eral parameters were steadly reproduced in all runs. These are 

(1) slightly undersolar nitrogen abundance, [N/C] = —0.1 ±0.1, 

(2) significantly lower ratios of Si/C and Al/C as compared to 
their solar values, [Si/C] = -0.3 ± 0.1, [Al/C] = -0.3 ± 0.1, 
and (3) extremely high overabundance of iron [Fe/H] > 1.5. 
The extremely high overabundance of iron as well as under- 
abundances of silicon and aluminium may be caused by non- 
equilibrium ionization. When an absorbing gas comes to equi- 
librium through cooling and recombining its ionization param- 
eter can be overstated due to different cooling and recombina- 
tion times of the observed ions stemming from density fluctu- 
ations (so-called "hot photoionization"). A higher density gas 
has shorter cooling and recombination times and hence comes 
faster to equilibrium. Thus, the sub-system A may be described 
as consisting of dense gas clumps that are already close to 
equilibrium (seen in low ionization absorptions) and ambient 
rarefied gas still far from equilibrium (responsible for most 
of C IV and N V absorptions). The equilibrium ionization pa- 
rameter should be lower, probably ranging between 0.001 and 
0.003. With this J7o, only a mild overabundance of iron, [Fe/Si] 
~ 0.2 — 0.3, would be enough to describe the observed intensity 
of the Fe ii lines. 



Thus, we estimate a metaUicity of > '^Zq in the higher 
ionized subsystem B and ~ 2Zq in the lower ionized sub- 
system A (both metallicities are referred to siUcon lines). This 
strong metalUcity difference is another argument in favor of 
non-equilibrium ionization in sub-system B. It probably has 
lower gas density than in A which leads to longer cooling and 
recombination times. We suggest that some time ago the whole 
absorbing complex {A + B) was either exposed to a much more 
intense radiation or shock-heated up to the temperatures when 
coUisional ionization becomes significant (Klein et al. 1994, 
2003; Levshakov et al., 2004). 
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Fig. 5. Same as Fig. 1 but for the Za,bs = 1.7103 system. The zero radial velocity is fixed at z = 1.71027. The corresponding 
physical parameters are hsted in Table 3, Cols. 3,4. The arrow in the Lya panel indicates pixels with non-zero intensities 



4.4. Absorption system at z^bs = 1.6838 

This system separated by ^ 12 450 km s"'^ from the QSO 
shows a partially saturated Lya and many metal lines with 
complex profiles (Fig. 6). Sill A1260, 1193. 1190 and Nv 
A1242 are blended with Lya forest absorptions, and at the po- 
sition of Si II A1526 a clear continuum window is seen. 

The apparent structure of the Lya profile suggests that 
this system can be divided into two subsystems: one at - 
250 km s-i < t; < -180 km s^^ (A) and another at -180 
km s~^ < V < 150 km s~^ (B) which were analyzed sep- 
arately. Calculations were carried out using both the ioniz- 
ing spectrum restored for the ^;abs= 1-7817 system and sev- 
eral power law spectra. Physical parameters hsted in Table 3, 
Cols. 5,6 correspond to the Zabs = 1.7817 background. The red 
wing of the synthetic profile of Lya in the subsystem B is cal- 
culated simultaneously with metal lines assuming a constant 



metallicity inside the absorber. The ionization parameter given 
in Table 3 for the subsystem A is determined from the observed 

intensity of C IV and the noise level at the expected position of 
the C II A1334 line and should be considered as a lower limit. 

The almost solar ratio of [Si/C] = 0.06±0.1 and a strong un- 
derabundance of nitrogen [N/C] < — 0.5 obtained for both sub- 
systems indicate that photoionization is probably close to equi- 
librium. However, slightly higher metallicity in the sub-system 
A and the remaining flux 0.04 ± 0.01 at the shallow bottom of 
Lya may indicate that hydrogen has not yet reached its equi- 
librium with the ionizing background (see sub-system B in the 
foregoing section). Although the absolute values of the abun- 
dances are uncertain (they depend on the assumed background 
and the mean ionization parameter), the solar to oversolar car- 
bon content, the ratio [Si/C] ~ and a significant underabun- 
dance of nitrogen are constantly reproduced in all trials. 
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Fig. 6. Same as Fig. 1 but for the 2:abs= 1-6838 system. The zero radial velocity is fixed at z = 1.68377. The corresponding 
physical parameters are listed in Table 3, Cols. 5,6 



4.5. Absorption systems at Zabs = 1-7365 and 
1.4978 

Other systems with unusually strong and complex C IV profiles 
are identified in the spectrum of HE 0141-3932. Unfortunately, 
they cannot be analyzed by the MCI because their hydrogen 
lines are unavailable, and we describe them only qualitatively. 
The system at Zabs = 1.7365 (separated by 7 000 km s^^ from 
the QSO) shows also a clear Nv doublet (Lya profile is 
blended). The apparent ratio C IV/N V is very much like that 
in the sub-system A at Zabs = 1.6838 and the Zabs = 1.7365 sys- 
tem probably has similar physical parameters. The system at 
z^t,g= 1.4978 {Av = 30 000 km s^^ ) exhibits a strong and 
complex Si IV absorption and a weak N V doublet as well as 
C I1AI334 (La is out of range), and in these features resembles 
the sub-system B at Zabs = 1-6838. 

4.6. The distance from ttie ligtit source 

The oversolar metallicity obtained in the absorbers described 
above place them in a class of associated systems. Furthermore, 
their extremely large radial velocities imply that they originate 
in gas ejected from the circumnuclear region of the QSO/host 
galaxy. This is also supported by the relatively strong emission 
flux in Fe II and Fe III lines and at the same time by the Fe II 
absorptions in the associated systems. Note that Fe II lines in 
absorbers with 7V(H l) < 10^^ cm~^ are extremely rare; prob- 
ably this is the first such detection. 



The distance between the absorbing cloud and the light 
source can be estimated from the photoionization model and 
the column density ratios of C II*/C II or Si Il*/Si II. For an ion 
in the interstellar (intergalactic) medium, the ratio of excited to 
ground-state population is equal to the ratio of the collisional 
excitation rate Qi^2 to the spontaneous transition probability 
A2^i (Bahcafl & Wolf 1968): 

The corresponding atomic data for C II* and Si II* are the fol- 
lowing: ^2^1 = 2.291 X 10^^ s^^, the excitation rate by col- 
Hsions with electrons at Tkin = 10^ K (^^2 ~ 1 x 10^'' cm'^ 
s-i for C II*, andAa^i = 2.17 x 10""* s-\ q\^^ ~ 2 x 10""^ 
cm^ s"^ for Sin* (Silva & Viegas 2002). Since colHsions 
with other particles have much lower excitation rates, we put 

We do not detect Si II* or C II* lines in the associated sys- 
tems, therefore clear 'continuum windows' at the expected po- 
sitions of C ll*A1335-7 and Si II*A1264.7 were used to set up- 
per limits on the column densities of C II* and Si II* (Table 3). 
For the 3 a upper limits on A^(Si II*) and A^(C II*), eq.Q pro- 
vides rio < 3 cm^'^in the 2;abs= 1-78, and 1.71 {A) systems, 
and ric < 15 cm^'^in the 2;abs= 1-68 (E) system- Since the 
degree of ionization in these systems is high (riH+ /f^H ^ 1), 
the upper limits on the total gas density nn are the same, i-C-, 
riH < 3 cm~'^ and nn < 15 cm~^ , respectively (the contribu- 
tion of the ionized helium is ignored since it has a small effect). 
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Table 3. Physical parameters of the Zabs= 1.7817, 1.7103 and 1.6838 metal absorbers toward HE 0141-3932 (zem= 1.80) 
derived by the MCI procedure (Umits are given at the 1 a level) 





Zato =1.7817 


^abs — 


1.7103 


^abs - 


1.6838 


Parameter 




subsystem A 


subsystem B 


subsystem A 


subsystem B 


(1/ 


(2) 


(3f 


(4)d 


(5) 


(6) 




S.IE-S" 


9.0E-3 


3.8E-2 


> 2.8E-2 


1.3E-2" 


ATh, cm 


4.2E17 


5.7E18 


8.6E17 


> 3.2E17 


2.5E18 


cTv. km s~^ 


20.8" 


55.0 


21.0 


19.2 


59.0 


o-y 


0.63" 


0.77 


0.62 


0.36 


0.52" 


"Zc 


3.4E-4" 


l.lE-3 


2.1E-3 


< 5.2E-^ 


2.9E-^ 


Zv. 


<1.2E-4 


3.1E-4 


3.2E-4 


< 5.0E-5 


2.8E-5'' 


Zms 


7.0E-5' 


1.4E-4 








Zai 


7.0E-6'' 


6.2E-6 








Zsi 


4.3E-5" 


7.4E-5 


3.0E-4 


< 8.5E-5 


4.8E-5" 


Zpe 


6.0E-5'' 


~8.0E-4 








[Zc] 


0.14 ± 0.08 


0.65 


0.94 


< 0.33 


0.08 ± 0.10 


[Zn] 


<0.15 


0.56 


0.58 


< - 0.3 


-0.48 ± 0.15 




0.30 ± 0.10 


0.6 








[Zai] 


0.30 ± 0.10 


0.3 








[Zs.] 


0.09 ± 0.08 


0.33 


0.94 


< 0.4 


0.14 ±0.10 


[Zfc] 


0.30 ± 0.15 


^^^1.5 








Af(Hl), cm~^ 


(2.6 ± 0.5)E15 


1.0E16'' 


(2.6 ± 0.5)E14 


(1.1 ± 0.2)E14 


2.5E15'' 


N{C II), cm~^ 


(1.64 ± 0.08)E13 


(1.5 ± 0.2)E14 


< 1.0E12 


< 1.6E11 


(7.3 ± 1.5)E12 


A''(C II*), cm~^ 


<5.0E11 


<4.3E12 






<1.5E12 


A'^(Mgll), cm~^ 


(1.20 ± 0.20)E11 


(5.3 ± 1.5)E12 








A'^(Si II), cm"^ 


(2.8 ± 0.3)E12 


(1.2 ± 0.1)E13 








iV(SiII ), cm 


<1.3E11 


<2.0E11 








iV(Feii), cm"^ 


(1.4±0.4)E11 


(5.5±2.0)E11 








7V(Aliii),cm-^ 


(3.6±1.0)E11 


(1.4±0.3)E12 








A''(Siiii), cm"^ 


7.5E12= 


(5.5±0.6)E13 


(2.3 ± 1.2)E11 


< 8.5E10 


(7.9±0.8)E12 


iV(Civ), cm-^ 


(2.2 ± 0.2)E13 


(1.7±0.2)E15 


(3.4±0.3)E14 


(3.4±0.2)E13 


(2.3 ± 0.1)E14 


Af(SilV),cm-^ 


(4.5±0.4)E12 


(7.7±0.8)E13 


(1.6±0.5)E12 


<3.0E11 


(1.2±0.1)E13 


Ar(Nv), cm"^ 


<1.4E12 


(1.8±0.2)E14 


(7.7 ± 0.8)E13 


(4.5 ± 0.5)E12 


(1.3±0.1)E13 


(T),K 


0.8E4 


1.2E4 


1.3E4 


1.8E4 


1.6E4 



"■''Uncertainties: " ~ 15 - 20%, and ~ 30%. 

°The value of N(Si ill) is calculated on base of other lines. 

'^Physical parameters (Uo . Nn ,a^.ay and metallicities) are only illustrative since estimated for the 

photoionization equilibrium whereas the absorption system is in non-equihbrium. 
''Estimated under the assumption of constant metaUicity inside the absorber. 

^Zx = iVx/iVn; [Zx] = \og{N^/Nn) - \og{N^/Nn)Q. 



To estimate the distance, the QSO continuum luminosity 
at the Lyman Umit Ci,^ must be known. Absolute spectropho- 
tometry is not available for HE 0141-3932. Therefore we esti- 
mated the intrinsic luminosity C^^ from {i) the comparison of 
the QSO B magnitude (assuming the QSO spectral energy dis- 
tribution is a MF-type, i.e., a — 0.5 in the range 912 A<C A <C 
1600 A) with the specific flux of a star having me = 0.0 
outside the Earth's atmosphere (4.4 x 10~^° erg s~^ cm~^ 
Hz~^) and from (ii) the empirical formula given by Tytler 
[1987, Eq.(17)] which is obtained from a fit to {f^, rriy) data 
on over 60 QSOs with a wide range of redshifts. The observed 
B magnitude of 16.09 (corrected for extinction) translates to 
the flux /^.(4400 A) = 1.6 x IQ-^e erg s'^ cm'^ Hz"! which 
in turn leads to the apparent luminosity near the Lyman limit 
- 8 X lO^i erg s"! Hz"! (ffo = 71 km s^^ Mpc^\ 
f2m = 0.3, and ^Ix = 0.7 are used to calculate the luminosity 
distance of 13 Gpc). The second method gives C^^ ~ 5 x 10^^ 



erg s~^ Hz~^ which shows that both values are in reasonable 
agreement. In the following we use the second one since it is 
based on the observational data. 

Given the upper limits on nn, the distance from the QSO 
to the absorbing cloud can be calculated from the estimated 
ionization parameter U which is defined as 

f/=j«ai^==*, (2) 

where 

Q(H0)= / -^dv (3) 

is the number of hydrogen ionizing photons emitted per unit 
time by the central source, c is the speed of Ught, Vc is the 
frequency of the Lyman continuum edge, and riph the corre- 
sponding density of ionizing photons. 
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With the estimated Lyman continuum luminosity Cy^ ^ 
5 X lO^i erg s"! Hz-\ one finds (3(H°) - 7.5 x lO^^^ photons 
s~^, assuming ~ C^^ {v/vc)~°', and a = 1 in the range 
V > i^c- A substitution of the numerical values in (|2ji provides 
ri.gg > 100 kpc, ri.7i > 280 kpc, and ri.78 > 450 kpc. 

5. Discussion and conclusions 

The characteristics of both the broad emission lines and asso- 
ciated absorptions of HE 0141-3932 can be best explained by 
the almost pole-on view of this quasar The weakness of emis- 
sion lines is then due to dilution by direct radiation from the 
accretion disk, whereas the velocity shifts of the associated sys- 
tems and their distances from the light source can be caused by 
the entrainment into the large-scale outflowing jet. However, 
powerful jets propagating to distances of ^ 0.5 Mpc imply the 
existence of significant radio radiation, but HE 0141-3932 is 
a radio-quiet QSO. An upper limit to its radio flux of 1 mJy 
translates into the intrinsic radio luminosity < 7 x 10'^^ 
erg s~^ Hz~^. This is a radio power of a typical FR I source. 
Jets in the FR 1-type structures are known to be subsonic or 
slightly oversonic with velocities 1 000-10000 km s^^ , heav- 
ier than and isobaric with the external medium (e.g., Hughes 
1991). These characteristics are in line with the suggestion that 
the observed absorptions originate in the entrained clouds. It is 
also known that FR I jets do not show an extended radio emis- 
sion. This means that radio observations of such jets, especially 
when they are seen at a small angle to the line of sight, need a 
high spatial resolution and high sensitivity. Nevertheless, the 
calculated distances of several hundreds of kiloparsecs seem 
to be 4-5 time overestimated since the typical jet length of the 
FR I object is below 100 kpc. We cannot explain the source of 
this discrepancy. It should be noted, however, that distances of 
hundreds of kiloparsecs for the associated systems are not ex- 
ceptional [see, e.g., Morris et al. (1986); Tripp, Lu, & Savage 
(1996); D'Odorico et al. (2004)]. 

The second fact that is hard to explain is the large veloc- 
ity excess between low- and high-ionization emission lines. 
According to a model of the quasar atmosphere (e.g., Elvis 
2004), low-ionization emission lines (Ha, O I, Fe II, Mg II, 
C III]) come from the outer region of the accretion disk which 
is well shielded from the central source radiation and optically 
thick for Lya. Lya and high-ionization emission lines (C IV, 
Si IV, N v) are formed in a cool phase of the wind arising from 
the inner parts of the accretion disk. The emission at 3400 A 
observed in HE 0141-3932 cannot be interpreted unambigu- 
ously, but if it is (even partially) due to Lya at z = 1.80, then 
its redshift and strength is inconsistent with this scenario since 
Civ and SilV are seen at z = 1.75. The same situation (i.e., 
ZLya > zciv) IS obscrved in PG 1407H-265. Apparently, a com- 
plex geometrical model of the broad emission line region is 
necessary to explain the observations. 

As already mentioned, the relativistic jet beamed toward 
us was discovered in the bright (B — 15.7) and radio-quiet 
PG 1407H-265 (Blundefl et al. 2003). There are indications 
that HE 0141-3932 may also have a similar small-scale rela- 
tivistic jet. The Lyman limit luminosity of HE 0141-3932 is 
-Ci/^ = 5 X 10'^^ erg s^^ Hz~^, whereas its radio luminosity is 



less than 7 x 10'^^ erg s^^ Hz^^. This implies either an unusu- 
ally flat radiation continuum (typically for QSOs is £^ oc i'^'^-^ 
at ly < 10^"^ Hz and oc i/^^-^ above), or that the apparent op- 
tical luminosity is Doppler-boosted due to the relativistic mo- 
tion of the light source. Another argument in favor of boost- 
ing comes from the metallicity estimations. HE 0141-3932 is a 
bright source. It is supposed that the QSO luminosity is deter- 
mined by the accretion rate which requires a large amount of 
circumnuclear gas. This, in turn, supposes a large mass of the 
contributing stellar population and, hence, a high metal enrich- 
ment of the accreting gas. For reference, all QSOs with lumi- 
nosities above 5 x 10"^^ erg s^^ Hz~^ in the sample of Dietrich 
et al. (2003) have metallicities Z > 3Zq with the mean value 
of 4-5 Zq. The metallicity of 4-5 Zq has been also measured 
in the associated system of a very bright QSO HE 0515-4414 
(Levshakov et al. 2003b). Gas in HE 0141-3932 has metal- 
Ucity Z w 1 — 2Zq, which supposes a relatively low stellar 
population involved in the enrichment, low accretion rate and, 
hence, low luminosity. Boosting can be caused by a relativistic 
jet seen at a small viewing angle. Taking into account the mea- 
sured equivalent widths of the emission lines, we may assume 
that the luminosity is amplified by ^ 10 times, i.e. the intrinsic 
luminosity of HE 0141-3932 might be only - 5 x 10^" 
erg s^^ Hz^^. 

Another issue of interest is the relative abundances obtained 
for the absorbing gas. The analyzed absorption systems show 
high iron content, [Fe/C] = 0.15 ± 0.1, [Fe/Mg] = 0.0 ± 0.1 
(zabs = 1-78), but at the same time nitrogen is strongly under- 
abundant, [N/C] < - 0.5 (zabs = 1.68). Although these values 
were estimated in different absorption systems they are rep- 
resentative for the bulk of circumnuclear gas for the follow- 
ing reason. The mass of the stellar population involved in the 
enrichment of a quasar's circumnuclear region is > IO^Mq 
(Baldwin et al. 2003a) and hence large metalhcity gradients 
and sharp discontinuities due to enrichment by only a few stars 
are unlikely. In the Hamann & Ferland (1993, 1999) models 
of QSO chemical evolution, solar metallicity is reached after 

> 0.2 Gyr and is characterized by a relative overabundance of 
nitrogen, [N/C] > 0, and an underabundance of iron, [Fe/C] 
< 0. Due to the delay of 1 Gyr in Fe enrichment expected from 
the longer evolution of SNe la which are the main sources of 
iron, the emission line ratios Fe Il/C IV and Fe Il/Mg II were 
proposed as a clock to constrain the QSO ages. However, in 
these models, large values of [Fe/C] and [Fe/Mg] are always as- 
sociated with a considerable overabundance of nitrogen, [N/C] 

> 0.3. We do not observe such behavior in our systems and, 
hence, cannot confirm this 'iron clock' . This is in line with the 
result of Matteucci & Recchi (2001), who showed that the time 
scale for enrichment by SNe la is not unique but a strong func- 
tion of the adopted stellar lifetimes, initial mass function and 
star formation rate and can vary by more than an order of mag- 
nitude. 
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Appendix A: Evaluation of thie spectrai shiape of 
thie iocai UV bacl<ground 

In general, the ionizing UV background is produced by the at- 
tenuated UV flux of QSOs and/or other sources (like young 
galaxies) and by the recombination radiation from intergalactic 
gas clouds. To estimate the shape of the UV continuum from 
the observed profiles of ions in the metal systems we use the 
response function methodology from the theory of experimen- 
tal design (see, e.g.. Box, Hunter and Hunter 1978, Chapt. 15). 

At first, the shape of the UV continuum in the range E > 1 
Ryd is to be parameterized by means of k variables (factors). 
For instance, the AGN-type spectrum can be described by a 
broken power-law defined by the following factors (see Fig. 4): 
/i - the power law index (slope) in the range 1 Ryd< E < f2 
Ryd; /2 - the coordinate of the first fracture; /a the slope 
between /2 and the point of the second fracture f^, /s the 
value of the break after (decimal logarithm of the inten- 
sities ratio) with the slope /g; /t the slope in the high UV 
range. Formally we fix the starting point of the X-ray break 
at £' = 100 Ryd and the slope after it at —1.5 since this spec- 
tral range affects very weakly the fractional ionizations of ions 
we are interested in. In this factor space, the Mathews-Ferland 
(MF) spectrum is represented by a point with the coordinates 
{-0.5, 0.24, -1.0, 0.61, -2.5, -3.0, -0.7}. 

At the start of the procedure we select a basic UV spectrum 
(e.g. power law, HM, or MF) and carry out the MCI calcu- 
lations using the fractional ionizations Ta,i' produced by this 
background. 

If the trial spectrum turns out to be inconsistent with the 
observed line intensities, we assign its parameters to a 'null 
point' in the factor space and vary the factors around this null 
point according to the chosen set of treatments^. Thus, for each 
treatment we obtain a new trial UV spectrum. 

The fitness of a particular spectrum is measured by the 
value of the response function TZ which is defined individu- 
ally for each absorption system in such a way as to ensure the 
self-consistent description of all absorption lines detected in the 
system. The best UV spectrum is that which provides a maxi- 
mum value of TZ within the constraints ^ 1 (P^r degree of 
freedom) for each individual line. The information needed to 
construct the response function is obtained in several test runs 
with the 'null spectrum'. 

The following example illustrates this step of the proce- 
dure. Let us assume that an absorption system exhibits lines 
of Si , Si IV and C II, C IV. Trial calculations with a basic UV 
continuum show that this continuum overpredicts the intensity 
of Si II and underpredicts Si IV whereas the behavior of carbon 
lines is opposite - underpredicted C II and overpredicted C IV. 

' If we define na,i as the total ion number density of element 'a' 
in the ith ionization stage (Wa — '^a.i), then Ta,i — nn^/un 
is the fractional ionization of ion {a,i}- Ror a given metallicity Z^, 

^ The treatment is a particular combination of values of the factors 
involved in an experiment. The set of treatments is called the 'exper- 
imental plan'. Experimental plans for any number of factors can be 
found in corresponding reference books (see, e.g., Johnson & Leone 
1977). 
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In this case it is conceivable to search for such a background 
that will produce a maximum value of the product Si iv/Si ii 
and C Il/C IV. The line intensities (and column densities) are 
proportional to the fractional ionizations and we can calculate 
the response function simply as 

n = ^x^. (A.1) 

Fractional ionizations for each trial UV background are taken 
from CLOUDY. Apart from the ionizing spectrum, further in- 
put parameters for CLOUDY are the ionization parameter U, 
metallicity and element ratios (e.g. [Si/C], [O/C]) which are 
estimated in test runs. 

After the response function has been evaluated at each point 
of the experimental plan, factor weights (influences) {a^jf^i 
can be estimated. For this purpose a standard polynomial model 
is used: 

k 

'R. = J^^aJi+(3, (A.2) 

where fi is the normalized and centered value of ith factor, 
fi = {fi — fo,i)/<^u - variability level of ith factor. Note 
that this form of TZ defines a response surface in fc-dimensional 
factor space. The coefficients and /3 are estimated from n 
(n > fc + 1) values of TZ computed at the plan points. The input 
of the nonlinear term can be deduced through the comparison 
of the estimate of [3 and the value of TZ^ at the 'null point', 
i.e. when = /o,^. The statistical significance of (7?,o — jJ) 
points to non-negligible factor interactions. Then additional 
points (treatments) should be added to the experimental plan 
to estimate the second-order terms. The optimal spectrum is 
obtained by moving the factor values in the direction normal to 
the response surface. For the absorption systems described in 
the present paper a linear model turned out to be sufficient. 

After the first iteration of the spectral shape adjustment is 
completed, we repeat the MCI calculations with the designed 
UV background and, if necessary, go back to the adjustment 
procedure until a satisfactory result (i.e. normalized ^ 1 for 
all profiles included in the analysis) is achieved. 



